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ABSTRACT 


A finite element based two dimensional transient 
heat transfer analysis has been performed to model the solidi- 
fication of steel and alundniura alloy castings. The work has 
been to arrive at optimal riser design criteria based on the 
solidification process modelling. 

The heat transfer formulation is based on heat condu- 
ction through the phases. The natural convective flow inside 
the liquid region has been neglected for the sake of simplicity. 
The heat release due to phase change in the solidifying mushy 
region is accounted through modified specific heat of the metal 
in the temperature range of solidification. 

Eight noded isoparametric elements have been used 
for discretizing the solution domain. Quadratic interpolation 
for temperature has been employed within each element. Gener- 
alized heat transfer boundary conditions of natural convective 
or radiative heat loss to the environment and heat conduction 
loss into the sand-base have been included. The transient 
temperature variation has been detained by marching in time 
using an implicit finite difference procedure. 

Based on the temperature results of the FEM analysis, 
isotherms were plotted within the casting at different times. 
From this the formation of isolated liquid regions which are 
cut-off from the riser are identified. If any such regions are 



predicted within the given casting configuration, the casting 
is termed unsound. Using this technique the optimality of 
riser dimensions and riser placement have been investigated 
for chunky and bar castings of steel and aluminium alloy* 

The results show good agreement with available classical 


works of riser design* 
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CHAPTER 1 


INTRODUCTION 


1 .1 GENERAL BACKGROUND 

Metal casting is an industry which involves high annual 
investment all over the world* An improvement in the quality 
of the castings and reducing the number of rejects would result 
in a large amount of saving for the casting industry. 

For manv years the production of sound castings was 
largely a matter of trial and error. In the early seventies, 
attention was focussed towards the scientific aspects of the 
production and design of castings. One such approach constituted 
the studies of the heat transfer mechanisms by which castings 
solidify and the rates at which solidification takes place. 

Of late, research of a more practical nature which is concerned 
with the design and placement of risers for castings has been 
carried out using modern computer simulation and graphical 
techniques. The present study deals with some of the practical 
aspects which need to be considered for the production of sound 
steel and aluminium alloy castings. 

The physical process involved in alloy solidification 
are primarily thermal in origin. A basic knowledge of the 
thermal phenomena will be helpful in meeting the various design 
requirements in a scientific way. The present study is aimed at 
developing an appropriate thermal model of the solidification of 
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sand-mould castings of alloy metals. It i s also intended to 
calculate the rate of solidification in the casting as a 
function of time and predict the occurrence of isolated liquid 
regions whose shrinkage due to solidification cannot be 
compensated by liquid fed from riser. By appropriately modify- 
ing the riser design, it is desired to eliminate casting defects 
such as shrinkage cavities with the help of the computer model. 


f *2 LITERATURE SURVEY 

Many investigators have contributed to the theoretical 
and experimental investigations on alloy solidification within 
a casting, over the past two decades. With regard to riser 
design, simplistic models exist which have led to the establish- 
ment of the current practices in the casting industry- However, 
very precise estimates of the riser dimensions and locations 
are not possible through these analyses. Limited work is also 
available on computer aided design of risers using modern simu- 
lation and graphical techniques. Here a brief survey of the 
available literature on alloy solidification and riser design 
is presented. 

Classical analytical solutions were derived by Neumann 
and Stefan for the propagation of a one-dimensional solidifi- 
cation front through a semi -infinite medium (1) considering only 
heat conduction within the solid and liquid phases. These 
solutions have been extended for a linear variation of thermal 
conductivity and for including radiative heat transfer in the 
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medium by (2,3) . Effects of buoyancy driven flov upon 

the rate of solidification have been considered by (4,5) 
for simple geometries, using the finite difference solution 
technique. These works have the limitation that they are 
not applicable for complex geometries and boundary condi- 
tions for which the solidification front may not be planar. 

Ostrach and Kroger (6) have considered the role of 
natural convection in the analyses of the steady state phase 
change associated with the continuous casting. While the work 
mentioned above concentrates on the effects of fluid motion 
an attempt to study the exact nature of the fluid motion 
itself has been made by Szekely and Stanek (7 ) * But the energy 
equation they considered did not contain any velocity terms. 
Hence the heat transfer process in not affected by buoyancy, 
conduction was the only transport phenomenon. 

Sparrow, Patankar and Ramadyani ( 8) carried out the 
analysis for two dimensional melting in a cylindrical enclosure 
taking natural convection into account. Their results were 
much different from those of the earlier conduction based 
solutions. Ramach andran, Gupta and Jaluria(9) studied the 
solidification in an enclosed region with natural convection 
in the melt for various heat transfer boundary conditions. 

They used an alternating direction implicit finite difference 
technique along with successive relaxation. 

Jaisuk Yoo and Rubin sk y ( 10 ) , have analysed the two 
dimensional solidification process with natural convection in 
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the melt, using FEM. Similar to Ram achandran, Gupta and 
Jaluria ( 9 ) , they considered a simple rectangular enclosure. 
They solved the governing equations using a novel "front 
tracking" Finite Element Method for predicting the interface 
movement . 

J.A. Sekhar, S. Kou and R. .Mehrabian (11) considered 
a model for the stationary transient heat flow problem during 
rapid melting and solidification of an alloy with a mushy 
Zone in which the heat of fusion is absorbed and liberated 
over a range of temperatures. A one-dimensional heat flow 
model which accounts for the heat of fusion was recently 
proposed by Hsu et al (12 ) . This model employs the approxi- 
mate integral technique coupled with the numerical technique 
developed by Murray and Landis (13) for the melting and solidi- 
fication of a pure metal* 

As regards the conventional practice of computing the 
optimal size of riser that will feed a given casting the methods 
adopted in the casting industry are largely based on the 
pioneering works of Chvorinov (14), Caine (15), Bishop (16) 
and Merchant (17). 

According to Chvorinov (14) the freezing time of a 

2 

casting is governed by its (V/A) ratio, where V is the 
volume and A is the surface area of the casting. Once 
(V/A)^ ratio of the casting is computed, the (V/A)^ ratio 
of the riser is selected such that it is a little larger 
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(about 2C^ ) than that of the casting so that the riser solidi- 
fies later. Caine (15) , who used Chvorinov’s relation, 

defined a new terra called relative freezing time, which is the 

o o 

ratio of (A/V) ratio of casting to (A/V) ratio of riser. 
Considering the shape variable, Bishop (16) restated the 
Chvorinov's rule and presented a curve of (L+W)/T vs. the ratio 
of the riser volume to the casting volume, which relates the 
length, width and thickness of a sound casting to its freezing 
ratio. Based on these works Merchant (17 ) proposed a general 
equation for riser size calculations. Although, these methods 
give reliable results they are time consuming. Large errors 
are likely to occur with complicated castings (since the effe- 
ctive surface area of castings which contain edges and corners 
is substantially different from the real area). Also detailed 
information on the shrinkage cavity size and its location are 
not available from solutions based on the above methods. 

However, in recent times, the F.E.M. simulation tech- 
nique has considerably simplified and improved the casting 
design procedures. Prasanna Kumar et al ( 1 8 ) have attempted 
F.E.M. formulations of 1-D freeze wave without assuming freeze 
wave equations, and a further extension of the model to a 2-D 
situation. The basic F.E.M. procedure to determine solidifi- 
cation times has been formulated using 2-D heat transfer 
conditions. 

In their papers, Imafuku and Chijiiwa (l9 ,20) have 
presented a method vrtiich predicts the shapes of shrinkage 
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cavities and porosities in steel castings with sand moulds. 

The method is analytically formulated in terms of the funct- 
ions of solid fraction and shrinkage cavity ratio, by means 
of mass conservation law. A close agreement between the 
results of the experiment and the calculations with regard 
to both temperature and shrinkage cavities have been demons- 
trated • 

The work done by Abis (21 ) describes the use of 
computer methods to obtain a good description of the solidi- 
fication pattern in an aluminium casting and has been tested 
by means of a numerical/experimental work on a 3-D mould- 
metal system. A FEM simulation adopting non-moving mesh was 
used to calculate the transient thermal field during the phase 
transition . 

Recently G.L. Shekar (22) has considered the interactive 
graphics approach coupled with the Finite Element simulation for 
predicting the occurrence of cavities inside an alloy casting 
for a two-dimensional situation. He demonstrated the validity 
of his technique by comparing with the known empirical results 
of the risering practice for castings of standard shapes. 

However the heat transfer boundary conditions in Shekar *s work 
were not generalized for the sake of simplicity. 

In the present work, an approach similar to that 
of Shekar has been considered with more realistic heat 
transfer boundary conditions on the external boundary of the 
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casting. Also, a computer software has been developed using 
8~noded isoparametric curvilinear elements which are capable 
of providing more accurate solutions for the temperature field 
and can handle complex geometries, while in the work of Shekar 
less accurate 3-noded triangular elements were used. Further 
the present approach is more suitable for the inclusion of 
buoyancy-driven flow effects from solidification point of 
view, 

1 *3 OBJ£CTI^/ES OF THE PRESENT WORK 

The present study deals with alloy solidification within 
a sand mould with aporopriate heat transfer boundary conditions. 
The main objective is to determine the riser size and location 
for producing sound castings thereby eliminating the costly 
trial and error method. 

The FEM solution technique is used to solve the 2-D 
transient heat transfer equation. From the predictions of the 
analysis, the solidification process within both the casting 
and the riser can be visualized by plotting the isotherms at 
different times in a graph paper. From these patterns the 
isolated liquid regions which cannot be fed for shrinkage 
compensation from the view can be observed. If isolated 
regions are identified then the size of shrinkage cavities 
can also be predicted. Once the occurrence of isolated 
regions is identified, it can be eliminated by changing the 
size or location of the riserfs) thereby achieving sound 
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castings. Thus, an optional riser design through detailed 
heat transfer analysis f or an alloy casting is the objective 
of the present work. 


^ *4 organisation of the press rfT WORK 

The mathematical modelling of the problem has been 
dealt in Chapter 2. The problem description, the assum- 
ptions made and the governing equations are discussed in 
this chapter. The FEM formulation has been dealt with in 
Chapter 3. The results and discussions have been made in 
Chapter 4 which summarises the stud/ conducted on both chunky 
and bar castings. Conclusions and suggestions for future work 
are covered in Chapter 5. 
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CHAPTER 2 


MATHEMATICAL MQDELLIi^ OF THE PROBLEM 


A clear understanding of the mechanism of solidifi- 
cation and cooling of liquid metals and alloys is essential 
for the production of sound castings. During alloy solidifi- 
cation, many important characteristics such as the crystal 
structure and the alloy composition at different parts of 
the casting are decided. Moreover, unless a proper care is 
taken defects l^ke shrinkage cavities will occur. Here we 
shall briefly discuss the salient features of the process of 
alloy solidification. 

2.1 MECHANISM OF ALLOY SOLIDIFICATION 

An alloy, unlike a pure metal, does not have a sharply 
defined freezing temperature. Instead, it takes place over a 
range of temperatures. During the process, the solid phases 
separating out at different temperatures possess varying 
compositions. Due to all these facts, the direction of crystal 
growth in an alloy depends on various factors, such as 

''i) the composition gradient within the casting, 

(ii) the variation of solidus temperature with composition, 

and 

the thermal gradients within the mould. 


(iii) 



Each of these factors are discussed by considering the 
example of a solid solution alloy whose phase diagram is 
shown in Fig . 2.1. 
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Let the liquid alloy have the composition (of B 
in a). Also let be the freezing point of pure metal A, 
and and respectively the liquidus and solidus temper- 

atures (the upper and lower limits for the solidification 
range) for the alloy of composition Cq* As the liquid alloy 
is cooled down to temperature solid starts to separate 

out. The concentration of B in these solids is only ( C^) 
as is evident from the Fig. 2.1. As a result, the concen- 
tration of B in the liquid, near the solid liquid interface, 
increases to a value more than C^* Figure 2.3 shows this for 
the situation where the solidification front has progressed up 
to some distant d from the mould face. 

Now, considering two points P and Q within the liquid 
alloy, P being just beyond the solid— liquid interface as 
indicated in Fig. 2.3. The solidus temperatures corresponding 
to the compositions at P and Q are Tp and Tq, respectively 
(see Fig. 2.1). Let Tp and Tq be the actual temperatures at 
the points P and Q. Tq is greater than Tp due to the thermal 
gradient within the casting. If both Tq and Tp lie in the 
range Tp to Tq, then the liquid at Q is supercooled, whereas 
that at P is not. This implies that crystallization starts at 
Q sooner than at P. If this difference is very prominent, 
the columnar growth of crystals (See Fig. 2.2) starting from 
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the mould surface is hampered which is undesirable. The 
crystal growth in such a situation appears as in (Fig. 2.4). 

Thus a dendritic crystal structure is produced. If the 
crystallization at Q gets completed before it starts at P 
(due to a very small thermal gradient) then randomly-oriented 
crystals may anpear inside the casting. Moreover, the presence 
of solid crystals ahead of the solid-liquid interface makes 
feeding of the liquid metal more difficult. This also implies 
greater risk of having voids within the casting. 

2.2 TEMPERATURE DISTRIBUTION OF DIFFERENT LAYERS WITHIN TH£ 

c.^tin6 

The heat rejected by the liquid metal is dissipated 
through the mould wall* The heat, released as a result of 
cooling and solidification of the liquid metal, passes through 
different layers. The temperature distribution in these layers, 
at any instant, is schematically shown in Fig. 2.5. The thermal 
resistances '4iich governs the entire solidification process are 
those of the liquid, mushy, solidified region, the metal- 
mould interface, the mold, and the ambient air. These six 
different reqions are indicated by numbers 1 to 6 in Fig. 2.5. 

A detailed heat transfer analysis is required to predict the 
variation of temperature through various layers and the formu- 
lation of such a thermal model is described below. 
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Fig 2.1 Equilibrium Phase Diagram of A-B Alloy. 
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Fig. 2.2 Development of Columnar Crystals. 
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Distance From Mould Face 

Fig. 2.3 Variation of Concentration of B with 
Distance from mould Face. 



Fig. 2.4 Dendritic Crystal Growth Structure. 
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2.3 THERMAL MODEL FOR ALLOY SOLIDIFICATION 

A transient 2-D heat flow model is developed for alloy 
solidification with a mushy zone in which the latent heat is 
absorbed over a wide range of temoeratures . The equations 
and solution methodoloqv developed are qeneral enough to be 
aoplicable to a variety of other solidification processing 
problems. A 2-D situation has been modelled for the sake of 
simplicity. 

The geometry of the sand mould, casting and the sand- 
base on which the heat transfer analysis has been performed is 
shown in Fig. 2.6. The boundaries of the solution domain are 
subjected to different modes of heat transfer conditions such 
as convective heat transfer, specified temperature and in some 
cases radiative heat transfer. As a result of these boundaries 
conditions and heat dissipation from the melt, temperature 
gradients arise within the melt and the sand mould which leads to 
solidification ©f the liquid. As the solidification progresses, 
three distinct regions are observed, namely: 

fi ) liquid 

( i i ) mu sh y 

(iii) solid 

These zones have different thermophysical properties. Typical 
values of the properties of these three phases for the steel 
and the aluminium alloy considered in the present work 
listed in Table 2.1. 


are 



1,9,10 - Constant Temperature (T amb.) 
2,3,4, 6, 7, 8 -Radiation/Convection 
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Initially the molten alloy is considered to be at a 
temperature greater than its melting temperature. This addi- 
tional heat above that necessary for melting is called super- 
heat, This superheat is necessary to provide the fluidity to 
the liquid vJhile pouring and It allows the melt to be trans- 
ferred into the mould freely. Thus the melt can access all 
the locations within the mould without starting to freeze even 
after coming into contact will the cold wall of the mould. 

The main assumptions that have beai invoked in the present 
investigation are: 

(a) The density, conductivity and heat capacity of the solid, 
mushy and liquid regions are independent of tamperatures 
but different from each other 

(b) For the mushy region an average thermal conductivity and 
a high heat capacity (due to latent heat) are used,^ 

(c) The effect of natural convection in the liquid zone is 
neglected as the rate of heat transfer is enhanced only 
by a small factor, 

(d) It is assumed that the sand mould is kept on a very large 
sand base. The sides and bottom of the sand base are 
assumed to be at roan temperature. The top surfaces 

of the sand bajse and the mould along with the sides of 
the mould are subjected to natural convective heat loss, 

4- These assumptions are not very realistic during alloy 

solidification since the solid and liquid fractions are not 
throughout the mushy region. Hov;ever, the major effects of 
property variation are still accounted for, by the consider- 
ation of average properties for each phase. 



(e) 


The melt poured into the mould has a prescribed 
degree of superheat* 


With these assumptions the 2-D transient heat transfer 
eguation governing the alloy solidification process is 
given by: 


P 

i 



3 t 


K. 




( 2 . 1 ) 


where the suffix 'i' denotes different phases such as liquid 
mushy, solid and mould* The above equation (^*1) is solved 
subject to the following boundary conditions: 


fi) On the surfaces 2, 3, 4, 5, 6, 7 and 8 (Refer Fig. 2 

where convective heat loss occurs 

(T - T-) = - for t ^ 0 (2.2a) 

C f o X 

(xi) For temperature specified boundaries 1, 9 and 10 

T = 7 , for t > 0 (2*2b) 

^ ^ arab ■" 

(xii) On surface 5, insulated boundary conditions is used 
since the top surface of the riser is sprinkled 
with some exothermic powder to prevent heat losses 
through it* Therefore, for this boundary 
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The case where both convective and radiative boundary 
conditions are specified on the surfaces 2, 3, 4, 6, 7 
and 8» 

(h_ + h^)(T - T-) = -K ~ for t > 0 (2.2d) 

re x on 

where 

hj. = o e ( I + T^) (T^ -j- T^) 

The initial conditions are: 

T = Tj_j+ at t = 0 within the casting and the 
riser (2 •26) 

T = at t = 0 in the sand ciculd and the 

sand base (2»2f ) 

For castings of various shapes placed inside send moulds t 
the heat transfer analysis is carried out numerically within 
the entire region of the casting, mould and sand-base shown in 
Fig* 2*6. After obtaining the temperature field as a function 
of time from such an analysis it is linked to the optimal riser 
design as described in the next section* 

2 *4 RISETi DESIGN AM D PL ACEMEN T 

VJhen risering a casting, it is necessary to design 
a riser with sufficient dimensions so that it will be the 
last portion to solidify and thereby contain enough liquid 
metal to suoply the shrinkage demand of the casting. Here, 
the conventional risering techniques are briefly discussed* 


(a) 

(b) 
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To properly determine the riser repuirement s of any 
castina the solidification characteristics of that casting 
must be evaluated. First, maximum feeding distances are 
calculated to determine riser locations. This can be done 
via the empirical rules (6) for the maximum feeding distance* 
Then the solidification time and the shrinkage volume of both 
the casting and the riser should be considered. 

"Tie solidification time of a given casting is commonly 
evaluated using the solidification modulus approach. The soli- 
dification modulus of a casting is the ratio of the casting’s 
volume to cooling surface area. Thus, 



According to Chvorinov's Rule ( 14 ) , 

t = (K (2.4) 

c 

where 

t = solidification time of the casting 
K = constant of proportionality (mould constant). 

Therefore, shapes having equal moduli have same solidification 
time when cast against the same molding medium. 

To produce a shrinkage free casting the riser is^ 
designed to have a larger modulus than the casting it is 
required to feed- This is to ensure that the riser solidifies 
after the casting. In practice. 
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= 1.2 (2.5) 

The 20?^ conspensation factor is necessary because 
of the nature of riser solidification. During the solidi- 
fication process the riser is constantly losing volume 
and increasing its top cooling surface area. Therefore 
the riser modulus decreases on solidification. Upon the 
completion of solidification the residual riser modulus 
should not be lower than the casting modulus in order to 
get a sound casting. 

The information on the amount of molten alloy needed 
from the riser is used only to compensate for the shrinkage 
that takes place from the pouring time till solidification 
is completed. Depending on the alloy, the percentage of 
this shrinkage varies from 2*5 to 7.5. Thus, the use of a 
large riser volume (to ensure large solidification time) is 
uneconomical. So, a riser should be designed with minimum 
possible volume while maintaining a cooling rate slower than 
that of the casting. 

The riser is designed to have a low surf ace -to -volume 
ratio than that of the section it is to feed. The low ratio 
will reduce heat dissipation and hence the riser will remain 
liquid than an equal mass with large surface. For this 
Durpose, a sphere is the ideal shape, but a sphere is difficult 
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to produce with a prooer positioning in relation to the 
casting. It is more common to use a cylindrical shane that 
has a diameter larger than the part to be fed and a heipht 
equal to about one to one-half times the diameter. 

In the present study the riser equations are obtained 
from Ref. (17) and it is used to get an approximate initial 
guess value of the riser diameter. 

The riser equation is expressed as follows* 

25( 4a -t.ll = 2-D - -.-1,»2.75„ (2.6) 

A s V. 7 ?- 


An optirmira value of q (= H/D) has to be chosen to calculate 
the value of D from Eqn. (2.6). The H/D ratio of the most 
compact cylinder which has a minimum surface area for a given 
volume can be determined in the following way? 


A 

V 


tiDH + 2 , and 

7 E?H or H = ^ 

^ TtD 


Hence , 

A 


% D 




4V ^ (^)d2 
U ^2. 


For A to be minimum, 

|A=o or-^+iiD = 0 or = ~ * 

3U 
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Again, 


or 


£V 

Tl 


D^H 


H = D 


O'" 


(2.7) 


This shows that the optimum ratio of H/D is equal to 1 for 
a cylindrical riser. Equation (2.6) is also dependent on 
the percentage volumetric shrinkage along with H/D ratio. 

Two cases have been considered in the present study corres- 
ponding to that of an aluminium-alloy and steel. 

(a) Initial Riser Size Calculations for an Aluminium Alloy 

Percent volumetric shrinkage of aluminium alloy = 6.6. 

(i) H/D =, 1 

Using the above values in Eqn. (2.6) we get 

A - 8.25 V - G. 08415 V A = 0 (2.8a) 

c c c c 

(ii) H/D = 1.33 

Putting the above values in Eqn. (2.6) gives 

- 7.8375 DT^V - 0,06311 V A = 0 (2«8b) 

c c c c 

(iii) H/D = 2 

Using the above values Eqn. (2.6) can be written as 

D^A - 7.425 DT^V - 0,042075 VA " 0 (2. Sc) 

c c c c 

(B) Initial Riser Size Calculations for Steel 

Percent volumetric shrinkage of steel = 

(i) H/D = 1 


3.0 
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Putting the above -'alues in Eqn. (2.6) gives 

- 3.75 D^V - 0,03825 V A =0 (2.8d) 

c ^ c c 

(ii) H/D = 1,33 

Putting the above values in Eqn, \2,6) gives 

- 1,3125 - 0.0286 V^A^ = 0 (2,8e) 

c c c c 

To check the adequacy of the riser size for steel 
castings Caine's relationship was also used. It is known 
that solidification time is proportional to the square of 
the ratio/surface area. Caine's relationship, however, is 
based on the assumption that the cooling rate is linearly 
proportional to the surface area/volume. A typical risering 
curve for steel is depicted in Fig* 2,7, Here the ordinate 
of a point in the curve shows the volume ratio and the 
abscissa the freezing ratio and the subscripts c and r 
refer to the casting and the riser, respectively. For a 
given casting-riser combination if the point in Fig, 2.7 
falls to the right of the curve, the size of the riser is 
adequate. The equation for a risering curve is of the form 



where 

a = freezing constant of the metal " 

b = the contraction ratio from liquid to solid 
= a constant depending on the different media 
around the riser and the casting. 


c 



:> 



D 

O 

> 



Freezing Ratio( A /V]^^(A/ V)r 


Fig. 2.7 Proper Combinations ot Volume and 
Freezng Ratbs for Steel. 
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Typical values of these constants for steel casting are given 
in Table ( 2 . 2)0 

So far our discussion has been restricted to the ade- 
quacy of the riser size from the points of view of shrinkage 
and cooling rate. Another important aspect of risering is 
to ensure that the available liquid metal in the riser can 
be fed to the desired locations within the casting. In fact 
the thermal cradient, within the casting, during the last 
stage of cooling is the most important factor which controls 
the occurrence of casting defects. The minimum allowable 
gradient depends on the shape and size of the cross-section* 
Generally one central riser is satisfactory if, 

^max “ ^ steel plates (2.9) 

where *t* is the thickness of the plate 

(ii) ^raax “ steel bars (2.10) 

where ’t ’ is the thiclcness of the bar in cm. 

The feeding distance should be measured from the edge of 
the riser. If the feeding distance is not proper centre line 
shrinkage will occur and hence the casting obtained will be 
unsound. 

The problem of riser design and placement under consi- 

% 

deration was solved with the help of FEM technique starting with 
an initial guess for the riser as given by the above calculations. 



27 


As the solidification proceeds with time it is checked 
whether any isolated liquid region is formed within the casting 
which will result in shrinkage cavities. In case such a region 
is found, the riser size is changed until it vanishes. Optimi- 
zation of riser dimensions is carried out along with a search 
for the optimum placement of the riser such that no shrinkage 
defect occurs. Aluminium alloy and steel castings having a 
simple geometry such as bar or cube have been considered for 
convenience, although complex shapes can also be handled by this 
model easily. 

2.5 OPTIMIZATION OF RISER SIZE AND DIMENSIM 

The use of a riser of small volume than that required 
to feed a casting leads to shrinkage cavity. The use of a 
riser with a very large volume (to ensure long solidification 
time) is also uneconomical from the cost point of view. Hence 
a riser should be desiqned with minimum possible volume while 
maintaining a cooling rate slower than that of the casting. 

So initially to start with it was assumed that 
H/D = 1 and the diameter of the base of the riser was taken to 
be half the diameter of the casting base. In all the cases 
shrinkage cavities were obtained. The next step was to 
increase the riser size and this process was continued until 
all the isolated liquid regions disappeared from the casting. 
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It was observed that the optimal riser diameter was less 
as comoared to the riser diameter calculated with the 
help of equation (2.6). Thus, overdesign of the riser can 
be avoided using the present computer model. 
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TABLE 2.1 : 

THERMAL PROPERTI 
STEEL, SAND *MOUL 

ES OF ALUMJNIUM ALLOY, 

•D AND SAND- BASE 



. St ate 

Thermal 

conductivity 

(V//m°C) 

Specific 

Heat 

(j/kg°C) 

Density 
(kg/m^ ) 

ALUf/IMIU?,* 

ALLOY 

SCUD 

193.53 

1 162.04 

2690 

(Al-3Cu-4,5Si) 

LIQUID 

91 .96 

1254.00 

2560 


MUSHY 

142.12 

5095.42 

2625 

STEEL 

SOLID 

31 .00 

802.56 

7500 


LIQUID 

41 .80 

836.00 

7000 


MUSHY 

36.34 

42636.00 

7250 

SAND MOULD 

- 

0.7524 

1065.90 

1580 

SAND-BASE 


0*3254 

794.96 

1500 


30 


TABLE 2.2 




AT. 


■S -0 


Freezing Constant 


(a) 


Contract! on 
Rati o 
(b) 


(c) 


0.10 


0.03 


1 .00 


TABLE 2.3 s CALCULATION CONDITIONS 


Alloy 

Liquidus 
Temp.( °C) 

Solidus 
Tempi °C) 

Heat 
Trans • 
Coef f i' 
cient 

Alumi- 

nium 

625 

525 

5-8 

Steel 

1525 

1515 

10-15 


Ambient 

Temo. 

(°c) 

Melt 

Temp. 

(°C) 

Time 
Step 
(Sec. ) 

Said-base 
Bottom - 
Temp* ( C 

30 

675 

20 

30 

30 

1700 

20 

30 
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CHAPTER 3 

FINITE ELEMENT ANALYSIS 


3.1 INTRODJCTION 

The solution orocedure that has been adopted for 
the problem under study is the Finite Element Method (FEM). 

Some of the most important advantages are the easy handling 
of complex geometrical shape of the solution domain and 
the general manner in which the boundary conditions are 
introduced. Rather than requiring every trial solution to 
satisfy the boundary condition the boundary conditions may be 
applied after obtaining the algebraic equations for the 
assemblage. Since the boundary conditions do not enter into 
equations for the individual elements, one can use the same 
field variable model for both internal and boundary elements. 
Moreover the field variables need not be changed, when the 
boundary conditions are changed. In addition FEM is also success- 
ful in representing complicated material property variations 
that are difficult to incorporate in other numerical methods. 

FEM also provides a systematic procedure for the derivation 
of the weighting function for problems with non-linearity and 
arbitrary domain. In view of the above mentioned generality 
and flexibility provided by FEM solution technique, it has 
been employed to solve the present problem. 
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DESCRIPTIg'I OF THE ?£M SCUlT'^ON P^ 0C£3UR£ -ALHRKIN'S 
APPROACH 

In FEM, the governing differential equations are 
converted into an equivalent set of integral equations which 
sre minimized in some sense over the whole solution domain. 
The solution domain is divided into several small elements of 
known shape and the solution variable is approximated within 
each of these elements by suitable approximating functions. 

The integral governing equations are then evaluated within 
each element using the approximating functions and all such 
elemental integrals are assembled to obtain the final matrix 
equation to be solved. The above steps can be summed up 
as f ollows : 

(a) Obtain the equivalent integral equations for 
describing the problem. The boundary conditions of the 
problem would be transformed into appropriate boundary inte- 
grals. 

(b) Divide the solution domain into several elemeats 
and place chosen number of nodes within each element. Obtain 
the approximating functions for describing the field variable, 
within each element in terms of nodal values. A mesh-gener- 
ation program may be used to obtain the finite element mesh 

Fig. 3*1 . 

(c) Each governing integral equation over the entire 
solution domain is written as the sum of the integral over all 






Fig. 3.1 Finite Element Mesh for Solution Domoin . 
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elements. 

(d) A typical element is isolated and the approxi- 
mating functions for describing the field variables are 
obtained within that element in terms of their nodal values. 

For this element, firstly, the elemental integrals are eva- 
luated leading to the formation of an elemental coefficient 
matrix and elemental right hand side vector. 'I'hen the boundary 
integrals (if any) are evaluated to obtain the boundary coeffi- 
cient matrix and the boundary right hand side vector. These 
elemental and boundary coefficient matrices are assembled for 
each element to obtain the final form of the elemental coeffi- 
cient matrix. A similar operation is done for the elemental 
right hand side vector also. 

(e) The elemental coefficient matrices are then 
assembled for all the elements, giving rise to a global coeffi- 
cient matrix. The elemental right hand side vectors are also 
assembled leading to the formation of the global right hand 
side vector. 

(f) The matrix equations obtained in the previous 
step are solved to find the unknown field variables at all 

nodes * 

In the present work, the Finite element approximation 
based on the Galerkin Method has been adopted. 

3 .3 METHOD OF WEIGHTED RESIDUALS 

Weighted residual methods are, in essence, numerical 
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techniques which can be used to solve a set of partial diff- 
erential equations* The first step in the application of 
the weighted residual procedure is to assume that the unknown 
variable, 0, can be aporoxi mated over the whole domain by 

0=2 w. 0. (3.1) 

i=1 ^ ^ 

where are the interpolation functions described in terras 
of the spatial coordinates (X, Y) and 0^ are the values of 0 
at certain chosen locations. 

Consider the equation t 

£ (0) =0 (3.2) 

where 0 is the exact solution of the differential operator 

* 

An approximate solution, (S , given by eq* (3*1), gives rise 
to an error or residue R, while substituted into eq. (3.2), 
such that 

R = £ (0*) 

In the weighted residual methods, R is minimised over 
the whole solution domain in the integral sense with a set of 
arbitrary ’weighting functions', , such that over the whole 
domain, D, 

/ W, R dD = 0 (3.3) 

D i 

^ 1,2,....no 


where i 
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It the nutTiber of unknown '/alues of 't is n, then by 
usinq n linear Indeoendent weiahtinq functions, enough 
equations can be generated to solve all the unknowns. The 
only limitation, at this stage, placed on , is that they 
roust be positive, single-valued aid finite® 

The Galerkin’s Method is a very special case of the 
method of weighted residuals, where the weighting functions, 
(W^), are taken to be the same as the interpolation functions 
vrtiich are more commonly known as shape functions. 

Therefore, for Galerkin's Method, the residue minimi- 
zation equations can be written as: 

/ N £(25*) d D = 0 

D ^ 


0 =-- 2 N. 0. (3.4) 

j=1 ^ ^ 

3.4 FBA SOLUTION PROCEDURE TO THE ALLCY -CASTING PROBLEM 

A major objective in a conduction analysis is to 
determine the temperature field in a medium resulting from 
conditions that are imposed on its boundaries. Oice this 
temperature distribution is known it is possible to optimize 
the riser size. ^ 

The 2-D heat conduction equation is of the form 
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P t-J = K T + Q (3.5) 

p 0 i, 

assuming constant thermal conductivity. The heat generation 
term has been neglected .Therefore Ec[n.(3.5) can be written as 

P C_ -K T = 0 (3.6) 

P C? t 

The equation (3.6) will be satisfied if temperature T is 
exact. Let us assume an approximate value of temperature 

* f \ 

say T and putting it in equation (3.6) gives 

P Cp - K T* = R (3.7) 

* 

where R is the residue due to the approximate solution T . 

The next step is to minimize the residue R over the whole 
domain D by Galerkin 's method. This is done by choosing 
some weighting function say N-^ which is also the shape 
function. This gives 

// N. R dx dy = [ 0 ] (3.3) 

D ^ 

Substituting the value of R from Eqn . (3.7) into £qn. (3.8) 

gives * 2 ♦ 

// N, ( P C - K V T ) dx dy = [0] 

D ^ 

(3.9) 
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^ S D % |t" K T* dx dy =[0 j 

(3.10) 

The unknown variable T can be expanded in terms of its nodal 
values and shape functions asi 

* n 

T = S N. T. (3.11 ) 

j=1 ^ 

where n is the number of nodes per element. 

Vfith the help of the weak formulation the order of 
the second term in Eqn. (3.10) is reduced 


-K 

D 


2 ♦ 

\ ^ T dx dy 


K // (v N. .V T* - V ,(N, 2 T*))dx dy 
D ^ ^ 

(3.12) 


According to the divergence theorem (for a 2-D case) 

(^ . F) dx dy = 6 F . n dl (3.13) 

D - ” 

The last term is equation (3.12) using divergence theorem 
becomes 


K Z 

-D 


(N^ Z T*) 


dx dy 


K 


^ \ Z T* . n dl 

B 


// V , (N V T*) dx dy 
D ^ 


K ^ \(Z T*. n) dl • 

B 

(3.14) 


or K 
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where B denotes boundary integral. 

_ * 

^ ST 

But V T • n = 5 “ which is nothing but the norEal 

“ cn 

derivative to the boundary 


-K 


K-i T* dxdy = K 


fl 

D 


Z J T 


3T 


dxdy- /H, (K ~~)dl 
1 c n 


(3.15) 


or 


9 * 3N. ^ * 

-K // N. V ^ T dx dy= K *- 

D D 


a X 


/ N. (K~-) di 
1 c n 



oT 

) dx dy 

cY ^ 


(3.16) 


The probleE 


conditions. 





has provisions for including all types cf boundary 

q for heat flux boundary condition 

-X- 

h_ (T - Tjf) for a convective boundary condition 

C X 

h (T - T^) for a radiative boundary condition 
r ^ 

where =cc (T + + T^) 


In the present work, only s-jT.talrical-shapsd castinc:s 
have been considered ersd therefore, it is sufficient to solve 
only half the problen. This is done by consider:' ng a dorain 
as shown in Fig. (3.1) with the prescription of an insulated 
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boundary conditions on the Diane of syrncietry* 

For a convective boundary condition, 5qn. (3.15) 
takes the form 


-K N. ^ ^ T* dx dy = K 

D ^ 


3 N. 




ax 3 y 



) dx dy 


+ / N. h (t* - T. ) dl 

1 t 


(3.17) 


where C denotes the boundary with convection. Substituting 
Eqn. (3.17j in Eqn . (3.10), we get the final equation of 
the entire domain as 



// 

-D 


3 N. 


3 N 3N. 3N. 

3 X 57^ 37^*'^=' 


+ hr / N, N, dl 1 fT. ] = h T- [ /N, dl 3 

‘-0^3 3 f ^ i 

(3.18) 


In the above equation, the sununation has been dropped 
by adopting the convention that summation is implied if a 
subscript is repeated. 

In the present study, 8-noded iso-oarametric elements 
have been used (See Fig. 3.2). The term isoparametric implies 
that both the element geometry and the vari ation in the vari- 
able across the element are described utilising the same type of 
polynomial or shape function. The global coordinates (x,y) of a 
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point inside an elerrient can thus be expressed in temis ox 
nodal coordinates and sh aoe functions as it has been done 
for the field variables- Therefore, 

X = 2 N. X (3.19) 

i=1 ^ ^ 

y = S N y (3.20) 

i=1 ^ ^ 

The shape functions for the 8-noded isoparametric elements 
(Refer ?ig» 3-3) in terms of local coordinates t and T , are 
given by the following relations. 

For corner nodes: 


Nj = j (1 + ) (1 + ’■'i'"' )(5i5 + 

(3.21 ) 

For raid-nodes: 

Ni = 1 (1 - §^)(1 i h ) for S ^ =0 (3.22) 

% = I (1 -T-)^)(1 )for = 0 (3.23) 

Substituting eguations (3.19) - (3.23) in equation (3.18), 
the final form of the assembled matrix equations are: 

[GMASS] [Tj] + [GC0N][T. ]= [ GQEN ] (3.24) 
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Fig. 3.2 Parabolic Eight Noded Element Showing Typical 
Node Numbering and Associated Discrete Varia- 
bles. For Programming Anticlockwise Sequence is Used . 



Fig. 3.3 Typical Shape Function Variation for an Eight 
Noded Isoparametric Element . 
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where 

[ GMASS ] is the global transient matrix, 

[ GCON ] is the global conductivity matrix, 

[ GGEN ] is the global right hand side vector. 

In terms of shape functions or their derivatives, these 
matrices and vector are expressed asi 


GMASS 

m 

= E 

e=1 

PC ff 

D 

(e ) 

‘*i 

N. dx 

j 

dy 


(3.25) 

GCON 

m 

= s 

e=1 

[ K 

(e) 

3N. 

3N. 

. 


2 M 

— j. 

3 y 

3 N. 

• iXi ) dx dy 

3 y 



+ h 

/ 

(e) 

h "3 

dl] 


(3.26) 

GGEN 

m 

= S 
e=1 


) ^'^i 

dl 



(3.27) 


The equation (3.1S) represents a set of simultaneous 
first-order ordinary differential equations of the variable 
T and its time derivatives. 

Hence the finite difference technique has been adopted 
to approximate the time derivatives. This leads to the final 
matrix equations of the form: 

[A] [ T ] = [C ] (3.28) 

There are three major schemes for handling the time 
derivatives in the finite-difference technique, namely! 
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(i) Iniplicit sch“n:-6 

(ii) Explicit scheme 

(iii) Crank-Nicholson scheme* 

Applying a generalised scheme which considers a fraction 
of the explicit scheme and a fraction of the implicit scheme, 
Eqn. (3.24) can be rewritten as, 

K-f 1 

{ f GMASS + e A t GCCN ] } [ J 
= {[GMASS - (1 - 9) At GCON]}[Tj]^ 

+ [ GGEN ] A t (3.29) 


The three schemes mentioned above can be obtained 
from Eqn. (3.29) by assuming different values of the fraction. 

For instance, 

0=0 for explicit scheme 
© = 0.5 for Crank-Nicholson scheme 

e = 1.0 for implicit scheme. 

Since the implicit scheme is the most stable one and 
has no restriction as far as the time step is concerned, it 
has been used in the present study. 

For implicit scheme. Equation (3.29) simplifies to 

the forms 

K+1 

! [GMASS] + at [GCONIJ [Tj] 

= [GMASS ][ Tj ] " +[GGEN]At (3-30) 
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Equation (3.30J can be solved to obtain the temperature 
vector at the (^K+1 )th time level, knowing the temperature 
vector at the Kth time level. 

3*5 MATRIX SQLUTiav' PROCEDURE : 

Here the Frontal solution method has been used to solve 
the assembled matrix equations. It is a well-known fact that 
the method adopted for solving the assembled matrix equation 
has a significant bearing on the computer storage requirement 
and execution time i^en the total number of unknown variables 
is very large. It becomes very time consuming and the core 
memory problem is very acute when a large number of unknowns 
are solved by the matrix inversion technique. The Frontal 
solution technique can handle quite a large number of variables 
without the problem of large core memory and the computer execu- 
tion time is also quite less* 

The Frontal solution technique used in the present 
work is based on the direct Gaussian elimination procedure for 
solving the symmetric matrices ^ere the leading diagonal is 
always used as a pivot. For unsymmetric matrices encountered 
in a wide range of engineering problems, the most suitable pivot 
is not necessarily on the leading diagonal and a Frontal solution 
routine exists for off diagonal pivoting (23) . But since this 
tends to be more time consuming, the method used in the present 
work uses only diagonal pivoting and incorporates many features 
of the Frontal method for solving symmetric matrices. 
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The overall Frontal solution techniaue covers the 
follovinq steosi 

i) Formation of the eleniental matrices. 

ii) Assembling into a small, temporary global matrix 
of chosen size. 

iii) Introduction of known variable boundary condition. 

iv) Reduction of the global matrix using Gaussian 
elimination procedure. 

v) Back substitution. 

The primary objective of the Frontal method is the 
elimination of variables as soon as possible after their 
introduction, via aopropriate equations, into the global matrix, 
^//hen the contributions from all the elerisents to a particular 
node point have been assembled, the corresoonding variables 
associated with that node can be eliminated. The complete 
matrix is, therefore, never assembled, since the reduced 
equations can be eliminated from the core and stored on disc. 
The equations held in core, with the corresponding nodes and 
variables, are termed as the FRONT and the number of unknown 
variables in the front is termed as the FRONT WIDTH. 

The Front width changes continuously, since, if the 
nodal matrix contribution have been fully summed the corres- 
ponding equation reduction based on the diagonal pivot can 

be executed. 
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^or a non-symmet ric global matrix, a oreassianed 
global matrix core area is filled from contributing elements, 
the largest diagonal entry in the pre-assigned core is then 
found and used as the pivot in a direct Gaussian eliminate. on 
process. As the maximum, predetermined number of equations 
are eliminated, the corresponding reduced equations are vnritten 
on to disc and more elements and corresponding equations are 
taken into core. 

The equations, nodes and variables currently in core 
are termed active, those assiqned to disc as deactivated and 
those yet to appear in core as inactive. This is shown diagram- 
atically in Fig. (3*4)* 
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o -Inactive Nodes 
9 - Active Nodes 
• -Deactivated Nodes 


VTT^TTm Front 

Assembled Elements 


Y/::r777\ Next Element for Assembly 


Fig. 3. 4 Definition of Front 
Front Width. 


and Element Numbering for Minimum 
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CHAPTER 4 

RESULTS AND DISCUSSION 'S 

Results or the temperature field have been obtained 
for various casting and riser configurat ions using the r&M 
technique. The isotherms have been plotted to identify the 
zones of isolated liquid regions which lead to the formation 
of shrinkage cavities. For presentation of results chunky 
castings and bar castings of aluminium alloy and steel have 
been considered. Cylindrical risers are used in all the cases 
as it is usually used in metal casting industries. The validity 
of the results have been tested by comparing the predictions 
of the present work with the risering curve, for the soundness 
of a casting of given riser dimensions. These comparisons 
for the riser design parameters are discussed towards the end 
of the chapter. 

SOLIDIFICATION ANALYSIS 

In Figure 4.1 the process of solidification of a 
bar casting of alundniura alloy has been shown. The casting 
size is 85 mm length and 20 mm thick and the riser dimension 
is 30 ram dia x 40 mm height. It i s observed that the solidifi- 
cation begins from the ends of the bar and slowly proceeds 
towards the riser. This shows that the thermal gradient is 
appropriate and leads to directional solidification towards the 
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source of feea metal vhich is of utmost importance from the 
DOint oT view of the soundness of the casting. In ?ig» 4.1(a) 
the isotherm of 625 C is plotted. It shows how after a small 
time mushy zone is formed within the casting. The riser 
remains liquid at this time which shows that the riser can 
feed liquid metal to the casting to make up for the shrinkage 
which occurs as the casting solidifies. Figure 4.1(b) shows 
the end of solidification. The casting has solidified comple- 
tely. The isotherm of 525°C is shown. The mushy zone is within 
the riser which is desirable. It is clear from the figure 
that the riser dimension chosen in this case would give a sound 
casting though this is not the optimum riser size. Riser place- 
ment is appropriate as the length of the bar is within the maximum 
feeding distance. 

Figure 4,2 and Figure 4.3 show the difference in the 
solidification rate for bar castings of steel. In the former 
case only convective heat transfer loss from boundary is consi- 
dered while in the latter heat loss due to both convection 
and radiation from the exposed surface of the mould is taken into 
account. It is observed that the solidification rate is higher 
duci to the enhanced heat transfer when radiation along with 
convection is taken into account. 

4.2 RISER DESIGN AND PLACEMENT 
4.2.1 CHOICE OF RISER SIZE 

Figures 4.4 to 4.6 show the nature of the solidification 
process when the riser dimensions are changed- In all these 
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cases the analysis has been performed on a cubical casting 
of side 57 mm of aluminium alloy. 

In Figure 4,4 the riser is of size SO mm x SO mm. 

It is observed that the solidification starts from the bottom 
of the casting and the riser is the last oart to solidify which 
is desired. Though a sound casting will be produced with this 
type of riser, an optimum value of riser size is desired to 
reduce the cost and extra metal. It i s shown later that it 
is possible to reduce this riser size without affecting the 
soundness of the casting. 

Figure 4,5 shows the effect of changing the riser size. 
Here, the riser size is 30 mm x 40 mm. It is observed that 
the solidification rate is such that the riser starts solidi- 
fying earlier than that of the casting leaving an isolated 
liquid region in the casting. Hence the riser cannot feed 
liquid metal to the casting after its own solidification. 

This leads to the fonuation of a shrinkage cavities and hence 
the riser dimension needs to be modified. 

At this stage the riser dimensions of Fig. 4.5 were 
increased in steps and for each of them the temperature field 
was calculated and plotted graphically to predict the isolated 
liquid regions. This procedure was continued until the casting 
was devoid of isolated liquid regions completely. Here the 
riser dimensions were noted and it was found out that the optimum 
riser size was 40 mm x 80 ram (Refer Fig. 4.6) for a 57 ram cubical 


casting 
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4.2.2 RISER PLACEMEMT AND QPTI.VAL RISER DESIGH PARA.ViET£RS 

Figures 4.7 to 4.9 show the importance of the riser 
pla:ernent and its effect on the nature of castings. In all 
these cases the analysis have been performed for dumbbell 
shaped castings (cubical parts of 60 mm each and central 
connection 20 mm x 40 mm) of aluminium alloy. 

Figure 4.7 shows the isolated liquid regions during 
the progress of solidification. Two risers are used, with 
one at each extreme end of the casting. The riser dimensions 
are 20 mm dia x 50 mm height. It is seen that after 75 secs, 
the mushy zone is formed within the riser though some parts of 
the casting are still in liquid form. So the riser is unsuit- 
able and hence the riser placement should be changed. 

Since the intention is to produce a casting devoid of 
voids or shrinkage cavities, the risers in Fig. 4.7 were both 
placed at the centre of the cubical portions as shown in 
Fig* 4.8. Though the size of the isolated liquid region 
is reduced as compared to the earlier case is reduced, still it is 
unsuitable for this type of casting. 

Now the riser position is kept constant and its size 
is increased in steps until no isolated liquid regions are 
found. Two risers each of 50 mm dia x 60 mm height, have been 
found to be optimum for the casting configuration considered. 

In the final design therefore the risers may be placed as shown 

in Fig. 4.9. 


i 
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The results for these cases have been coTiOared with 
the findings of earlier works (24) and there is an excellent 
agreement between them. 

4*3 FEEDING DISTANCE 

Ati important aspect of risering is to ensure that the 
available metal in the riser can be fed to the desired locations 
v7ithin the castings. 

To check the accuracy of the predictions from the 
present model, comparisons were made with known results ( 24,25 ) 
on the maximum feeding distance of a riser. 

The general equations for maximum feeding distance 

(v ) has been given by Eqn. (2.9) and Eqn. (2.10). 

'*raax 7 -I 

For this aspect of the study only steel bar castings 
were analysed. Two cases were taken up* 

The first case considered is the casting geometry 
of Fig. 4.2 which shows a steel bar casting with thickness 
(t) 20 ram and length 90 ram, with the riser dimensions of 
30 mm dia x 40 ram height. The maximum feeding distance 
calculated from Eqn. (2.10) is more than half length of the bar 
provided in this present case. It is observed that the 
solidification proceeds from the extreme ends of the bar and 
no isolated regions within the casting. Hence tlje 
riser size, placement and feeding distance are appropriate. 

In the second case, i or the casting mentioned above 
the riser dimensions were maintained same along with the 
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casting thickness, but the casting length was increased to 
225 rnm as shown in !“ig« 4.10 so that half length of the bar 
exceeds the maximum feeding distance calculated from cqn*2»10. 

In the numerical results, isolated licuid regions are observed 
as expected. Hence the feeding distance should be brought 
within required limits to prevent the formation of shrinkage 
cavities. 

The results predicted in these cases have been well 
supported by the available results in literature (24,25) • 

Thus the present analysis provides a very powerful tool for 
the optimal design of risers. 

The validity of the predictions from the present work 
upon the soundness of castings are compared with those of riser- 
ing curve in Fig. 4.11. For the situations considered 

in Figs. 4.2 (Case l) and 4.10 (Case II). The volume 
ratio (V_/V ) and freezing ratio iA/V) /(A/V) have been 
calculated for these casting configurations and the points 
corresponding to these values have been placed in relation 
to the risering curve. It is observed that the point for case I 
falls in the zone of sound castings, while that for case II does 
not. These observations are in good agreement with the predi- 


ctions of the present work 
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4.2 Solidification of a Bar Casting of 0.15% Carbon Steel for Convective 
Heat Loss from the Boundary. No Centreline Shrinkage Area when Bar of 
Length Twice the Maximum Feeding Distance of Riser. 
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Fig, 4.6 Risering of a Cube Casting of Aluminium Alloy (92.5 % Al, 
3% Cu, 4.5% Ni ). Modified Optimum Riser Used. 
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Fig. 4.9 Risering of a Dumbbell Shaped Casting of Aluminium Alloy(92.5®/oAl, 
3% Cu,4.5®/oNi) Two Proportionate Risers One Each at the Centre 
of Each Casting. No Isolated Liquid Region Seen. 
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CONCLUSIONS AND SUlGESTICfIS 

5.1 CONCLUSIONS 

(l) For the particular type of castings analysed, 
changes in the riser dimensions and riser placement have a 
large effect in the soundness of the casting. 

^2) The present riser design procedure can give 
significant improvements in true casting yield. 

5.2 SUGGESTIONS FOR FUTURE WORK 

(1) The present problem can be easily extenaed to more 
practical situations of three dimensional geometry castings. 

(2) The phenomenon of natural convection v-ithin the 
liquid metal can be looked into and the enhancement of heat 
transfer due to this phenomenon can be precisely modelled. 

(3) Detailed variations in thermal properties can 
be considered and the analysis can be performed for various 

alloy metals. 

(4) Mass transfer aspects of alloy solidification 
and detailed flow considerations within the mushy zone 
through a porous medium approach, can be included. 
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